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Abstract: We compute the matrix elements of the energy-momentmn tensor between 
glueball states and the vacuum in SU(3) lattice gauge theory and extrapolate them to 
the continuum. These matrix elements may play an important phenomenological role in 
identifying glue-rich mesons. Based on a relation derived long ago by the ITEP group for 
J/ if) radiative decays, the scalar matrix element leads to a branching ratio for the glueball 
that is at least three times larger than the experimentally observed branching ratio for the 
/o mesons above IGeV. This suggests that the glueball component must be diluted quite 
strongly among the known scalar mesons. Finally we review the current best continuum 
determination of the scalar and tensor glueball masses, the deconfining temperature, the 
string tension and the Lambda parameter, all in units of the Sommer reference scale, using 
calculations based on the Wilson action. 
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1. Introduction 



SU(N) gauge theories in 3+1 dimensions have been studied by lattice Monte-Carlo tech- 
niques a long time [1]. Many low-energy dimensionless quantities are now known at the 
percent level for = 3, as we shall review. However, low-energy matrix elements of local 
operators have not been given much attention. This is mainly because the lowest possi- 
ble dimension for a local gauge invariant operator is four. For Monte-Carlo simulations, 
asymptotic freedom then implies that the amount of statistics has to grow like to main- 
tain fixed relative errors on the matrix elements [2] (a is the lattice spacing). This comes 
on top of the unavoidable cost of simulating a four-volume fixed in physical units. In 
particular it is very expensive to take the continuum limit of renormalized matrix elements 
of the energy-momentum tensor, as is well known to thermodynamics practitioners. In 
this paper we compute the matrix elements of the energy-momentum tensor between the 
vacuum on the left and a scalar or tensor glueball on the right. By using a locally (i.e. on 
the scale of one lattice spacing) smoothened gauge field, we are able to reduce the prefactor 
of the ^ cost function to a manageable size. 

While these matrix elements have been computed previously [3], our technology differs. 
In particular our lattice spacings are significantly smaller and we improve on the statistical 
accuracy of the matrix elements in the continuum by a factor of about two. This im- 
provement is due mainly to the reduced uncertainty on the non-perturbative normalization 
factors of the energy-momentum tensor. 

One of our motivations is the possibility to confront model predictions with the lattice 
data. Indeed there are QCD sum rule predictions for these matrix elements [4]. Fur- 
thermore, models of QCD based on the AdS/CFT correspondence readily predict many 
stables glueballs [5], but the spectrum of the SU(3) theory contains only two scalar and 
two tensor strictly stable glueballs, and this is likely true at any finite number of colors N. 
By contrast, the detailed properties of the low-lying states, beyond their mass, represent 
opportunities for unambiguous comparisons. The ability of the energy- momentum tensor 
to annihilate a glueball is a quantity of this type. Such matrix elements are somewhat 
analogous to F^, which determines the width for the vr'^ fi^v decay. 

Two other applications will be discussed in section 4. A long time ago the ITEP group 
derived an approximate expression [4] for the branching ratio for a J/ijj to decay into a 
scalar glueball in terms of the matrix element that we compute in this paper. We can thus 
roughly estimate the expected production rate and compare to the experimental branching 
ratios for scalar mesons. 

Secondly we show how the tensor and scalar matrix elements can be used to constrain 
the thermal spectral functions that determine respectively the shear and bulk viscosity of 
the plasma of gluons [6, 7]. 

In section 2 we define the matrix elements to be computed, and give the relation 
between the lattice observables and the continuum, relativistically covariant quantities. 
Section 3 describes the lattice calculation. In section 4 we compare our results with those 
of [3] and present the aforementioned applications. We end with a summary of the current 
knowledge of the low-energy properties of SU(3) gauge theory. 
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2. Definitions 



In this section we fix our notation, define the relevant matrix elements and show how they 
can be obtained from Euclidean correlation functions. 

2.1 Glueball matrix elements in the continuum 

Decomposing the energy-momentum tensor T^y into a traceless part O^y and a scalar part 

via T^y = 6^^ + jSi^uO, the explicit Euclidean expressions are 

eix) ^ (5{g)/(2g) F^,{x)F^,{x) e,y{x) ^ l5,yF^,F^, - F^F^^. (2.1) 

The beta-function is defined by qdg/dq = (3{g) = —g^{bo+big^+. . .) and bo = llAr/(3(47r)^), 
bi = 34iVV(3(47r)^) in the S\J{N) pure gauge theory. The gauge action reads = jF^^F^^ 
in this notation. 

We take over the notation of [3] and define the matrix elements 

{n\e{x)\S,p)r = s e-'P-"" (2.2) 
{n\eu{x)\T,p,a-)r = {n\'^{en - 022){x)\T,p,a+)r = t e-'P-'', p={0,0,P3), (2.3) 

where l^) and |T) respectively refer to the lightest scalar and tensor glueball states. The 
labels a± refer to the superpositions of helicity states \a+) oc | -|- 2) -|- | — 2) and |cr_) oc 

1 -|- 2) — I — 2) . We define s' and t' in the same way for the first excited glueball in each of these 
channels. The subscript 'r' (relativistic) indicates that the state has a Lorentz-invariant 
normalization, 

,{Upa\Uqa\ = 2p\2Trf6{p-ci)6aa', U = S,T. (2.4) 

2.2 Glueball matrix elements from Euclidean correlation functions 

We now discuss the extraction of glueball matrix elements from correlation functions in 
the Euclidean theory, set up in a finite (but large) spatial volume. Separately for the scalar 
and the tensor channels, we consider the correlation matrix 

A,,(r) = L-3 (0,(0) 0, (r))e, i, j = 0, 1, . . . iV^. (2.5) 

The subscript 'c' indicates that we are dealing with the connected part and r is the Eu- 
clidean time variable. Let us consider first the scalar channel. The operator with label 
is the definite-momentum projected local current, 

C)o(r) = y"d^x e^P-^^(r,x). (2.6) 

Operators 1 through Ng are definite-momentum glueball operators, typically extended and 
designed to have large overlaps on the lightest two states. In finite volume the spectral 
representation of (2.5) reads 

oo 

A,(r) = J2 (2.7) 

n=l 
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where states are normalized as in quantum mechanics, {n\ TTi^ — and have momentum 
p by momentum conservation. In the infinite volume limit, the connection between the 
first state |1) with the one-particle state \S,p)i introduced above is 

^2EpL^ |1) ^ \S,p),, (2.8) 
where = Mj + P^- In particular, at large Euclidean time r, 

^oo(r) = Fl e-^p- + F|, e'^p^ + O(e-^p-). (2.9) 

with 



Fs{L) = L-^/^\{0\Oo\l)\ and lim F5(L) = 

L— >oo 




(2.10) 



and similarly for the excited state |2) which in the infinite volume limit has energy E'^ = 
M|, + p2. 

In the tensor channel, similar equations apply as long as p is collinear with the po- 
larization axis. On the lattice, the equality between the two forms of Eq. 2.3 is violated 
by O(a^) discretization errors as well as finite-size effects. Here we use the second form 
exclusively, but we perform checks for both sources of systematic error. 

In the following we will only use p = 0. 

3. Lattice calculation 

We simulate the SU(3) gauge theory using the Wilson plaquette action [8] at three values 
of the bare coupling, (3 = G/Qq = 6.0, 6.2 and 6.408. This corresponds to values of the 
Sommer parameter ro/o = 5.368(35), 7.383(55) and 9.845(85) [9]. The lattice size is 
respectively 16^ x 24, 20^ and 28^. On the coarsest lattice, we also check for finite- volume 
effects by performing an additional simulation on a 20^ x 24 lattice. The glueball spectrum 
was determined rather accurately in [10] at these lattice spacings. At the smallest lattice 
spacing, we apply a conversion factor (given by the ratio of tq values) to convert the 
spectrum from /? = 6.4 to /? = 6.408. We use the standard combination of heatbath and 
overrelaxation [11, 12, 13, 14] sweeps for the update in a ratio increasing from 3 to 5 as 
the lattice spacing is decreased. The overall number of sweeps between measurements was 
also increased, from 8 to 12. 

We use the 'HYP-plaquette' discretization of the energy-momentum tensor, as de- 
scribed in [15]. We determine its normalization non-perturbatively, following the same 
strategy as in [15]. The normalization of the discretization based on the bare plaquette is 
fixed by lattice sum rules, as is well-known in the context of thermodynamics [16]. It is 
therefore sufficient to determine the normalization of the discretization based on the HYP- 
plaquette, which we employ here, relative to the bare plaquette. A straightforward way to 
do this is to match the pressure and energy density computed with either discretization at 
a specific temperature. The choice T = 1.21Tc made in [15, 2] insures that a large signal 
is obtained for {9) and that T~^(^oo) is already more than half its Stefan-Boltzmann limit. 
Secondly, this choice implies A^^- > 6 for /? > 6.0, so that large cutoff effects are avoided. 
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Xsigo) xigo) 

6.000 0.9951(77) 0.5489(68) 

6.093 0.9778(89) 0.546(14) 

6.180 0.976(12) 0.596(20) 

6.295 0.953(20) 0.563(28) 

6.408 0.985(42) 0.612(49) 





aMs 


aM*g 


aMT 


aM* 


6.000 0.7005(47) 
6.100 0.6021(85) 
6.200 0.5197(51) 
6.400 0.3960(93) 


1.167(25) 
1.038(15) 
0.929(10) 
0.690(18) 


1.0596(64) 
0.916(11) 
0.7784(79) 
0.5758(32) 


1.433(14) 
1.180(34) 
1.032(20) 
0.795(28) 



Table 1: Left: normalization factors for the HYP-smcarcd plaquettes relative to the bare plaquette. 
Right: the glueball masses from [10] relevant to this work. 



These relative normalization factors, denoted by Xs{go) and x{9o) respectively for 9 and 
^00) are given in table 1. The data is conveniently parametrized as (6 < 6/50 < 6.408) 

Xs(5o) = 0.9731 + 0.67(5^ - 6/6.18) (3.1) 
x(5o) = 0.5701 - 0.77(5^ - 6/6.18). (3.2) 

In both cases, the absolute error increases from 0.007 at /3 = 6 to 0.020 at /3 = 6.408. 
We extract dgQ^/dloga from the parametrization of log(ro/a) given in [9] and use the 
parametrization of Zi^go) given in [6]. 

We employ linear combinations of Wilson loops that project in the Af^ and E'^^ 
irreducible representations of the cubic group. They arc constructed from spatial links 
variables, using smearing and blocking as described in [17]. 

3.1 Extraction of the glueball matrix elements 

The glueball matrix elements can be extracted at sufficiently large r using the fit ansatz 
Mr)= Yl c«c(f'^e-^"- (3.3) 

n=l,2 

^ Fu = jcf^l , Fu' = jcf I, U = S,T. (3.4) 
A few remarks are in order: 

• the correlator ^00 is not included in the fit. 

• in our discretization the operator Oq is defined at half- integer times (in lattice units) ; 
thus for i > l,j = 0, rja takes half integer values in (3.3), whereas for i,j > 1 it 
takes integer values. 

• the glueball spectrum is already accurately known at the same simulations parameters 
from ([10], Table 7.1). We reproduce the values M„ in table 1 and treat them as a 
'prior'. 

Note however that since this information is exact within its quoted error, it is unnecessary 
to invoke Bayesian arguments to make use of it. Therefore we determine the fit parameters 
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by minimizing = Xi + Xii) with 

= E E E [^^^•(^) - ^*^(^)] C'-\.,r),(fc;r') [^feKrO - Akiir')] (3.5) 

i<i T=«ij T'=Ukl 

k<l 

^.^^(M^^ (3.6) 

n=l,2 ^" 

We use the full correlation matrix C in this fit. Autocorrelation effects in Monte-Carlo 
time were found to be negligible, so that C(^ijT-'j^(^i^iT-i-j could be taken to be the usual estimator 
of ^mlt[iAiji'T)Akli'r')) - {Aij{T)){Aki{T'))], where (.) are path integral averages and Ny^st 
is the number of measurements. We estimated the relative error on the eigenvalues of C 
and found them to be small, due to the high statistics of the calculation. This makes the 
inversion of C a stable procedure. 

We used two non-local operators in the fits, i.e. the indices i,j,k,l < 2 in Eq. 3.5. 
These operators were linear combinations of non-local operators in the relevant lattice 
irreducible representation designed to have a good projection on the lightest two states. 
Increasing the number of operators in the fit did not seem to improve the determination of 
the physical parameters. The results and minimized x^ of these fits are given in Table 2. 
The x^ arc of order unity, with a tendency to be slightly larger that one. This may be due 
(partly) to the fact that the standard estimator for the x^ we are using is an upward-biased 
one [18]. Beyond the value of the x^ it is important to check that the fit is stable against 
variations in the fit ranges Uij < r < Vij, and that no strong trend is seen in (r)/ (r) 
as a function of r within the fit range. Figure 1 displays this ratio in the scalar and in the 
tensor sectors for the simulation at the smallest lattice spacing. We see that the correlators 
^01 cind All have small deviations from the fit even outside the fit range. This gives us 
confidence that excited state (n > 3) contaminations are negligible in the fit range. This 
property is much less satisfied for the A02 and A22 correlators. Although the fit range 
for the correlators starts further out in r, this implies that the control over excited state 
contributions is much less good for Fs* and Ft* . We will therefore not include these matrix 
elements in our final list of results. 

We may use the effective mass and the effective matrix element 

ameff(t) = log (3.7) 
Aii[T + a) 

FuM-r - = ^01 - \^^^y/2a ' U = S,T, r = 2a, 3a, . . .(3.8) 

as a way to visualize the contributions of excited states to the physical quantities we extract 
(see Fig. 2). The effective quantities computed from the fit to the correlators (Eq. 3.5,3.6) 
are also shown as curves on the plot. The fit looks convincing. We find that, in general, 
the method of determining the mass and matrix elements from the effective quantities is 
less stable than fitting the correlators. A possible reason for this is that the histograms of 
ratios as in Eq. 3.8 can become arbitrarily non-Gaussian if the denominator is noisy (e.g. 
at large r), even if the correlators were perfectly Gaussian distributed. 
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Data / Fit ( Scalar, p = 6.408 x = 1 -79 ) 




Data / Fit ( Tensor, p = 6.408 x = 1 -34 ) 




Figure 1: The ratio of the scalar (top) and tensor (bottom) correlators to the fit, Eq. 3.3. The 
arrows indicate for each correlator where the fit starts. 



Finally, increasing the volume from 16^ to 20^ at /3 = 6.0 does not affect the matrix 
elements in a statistically significant way (see Tab. 2). An effect at the subpercent level is 
seen on some of the glueball masses, because they are so accurately determined, but this 
does not affect the discussion in the rest of this paper. 
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Effective Mass and Matrix Element 

5.5 I — ' — I — I — I — I — I — I — I — I — I — I — I — I — ' — I — I — I — I — I — I — I — I — I — I — I — ' — ' — |— 




2 I ■ ■ ■ I ■ ■ ■ I ■ 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



Figure 2: The effective mass and matrix element in the scalar sector at /3 = 6, 16^ x 24 lattice. 
The fit Eq. 3.3 is shown. The fit starts at = a ~ 0.19ro for the diagonal correlator and at 
a;o = |a Ri 0.47ro for the correlator between the local and the smeared operator. 





aMs 


aMs* 


aMr 




6.000 (20^) 

6.000 (le^) 

6.200 
6.408 


0.7032(16) 
0.6966(16) 
0.51918(96) 
0.39603(63) 


1.1733(37) 
1.1747(58) 
0.93010(30) 
0.6946(11) 


1.0587(15) 
1.0592(22) 
0.7783(20) 
0.5937(36) 


1.43261(53) 
1.4358(10) 
1.0399(22) 
0.8322(92) 



6.000 (20^) 5.05(06)(09) 6.53(30)(12) 1.33 3.63(26)(08) 8.49(78)(27) 1.29 

6.000 (16^) 4.95(03)(09) 6.21(19)(12) 1.36 3.41(11)(08) 7.72(40)(25) 0.94 

6.200 4.43(25)(09) 10.25(90)(22) 1.50 3.07(32)(09) 10.37(68)(29) 1.29 

6.408 4.55(36)(14) 9.75(54)(31) 1.79 2.30(33)(10) 11.04(38)(48) 1.34 

Table 2: The scalar (S) and tensor (T) glueball masses and matrix elements extracted from the 
fits. In the lower table, the first error is the uncertainty coming from the bare matrix element in 
lattice units, the second is the cumulated error of all the other factors entering the renormalization 
group invariant quantity. 



3.2 Continuum extrapolation 

Having obtained the matrix elements at three different lattice spacings, we can attempt a 
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Continuum Extrapolation of Glueball Matrix Elements 




(a / ro)-^ 

Figure 3: Continuum extrapolation. 



0.04 



continuum extrapolation linear in (a/ro)^, as illustrated on Fig. 3. We find 

Fsr^^^ = 4.13(40) , xVd-o.f. = 0.46 (3.9) 
FtvI^^ = 2.10(41) , xVd-o.f. = 0.99 (3.10) 

The extrapolations are satisfactory: the is small, the slopes of the extrapolations are 
small to moderate. In particTilar the continuum value is statistically compatible with the 
value from the smallest lattice spacing. 

We do not attempt an extrapolation for the first excited states. Indeed it is clear from 
Fig. 3 that the /3 = 6.0 data does not lie in the scaling region. The values from the 
lighter two lattice spacings are consistent with eachother, suggesting that cutoff effects are 
small beyond /3 = 6.2. Our best estimate of Fs/Tq^^ and Ft'Tq^^ are thus the values at 
the smallest lattice spacing, but one should keep in mind that our control over excited 
state contamination is far less good for these matrix elements; the induced systematic 
uncertainty is comparable or larger than the statistical error. 



4. Applications and summary 

In the following, we compare our results to previously obtained ones. Then we present two 
applications of the computed matrix elements. We end with a summary of our knowledge 
of non-perturbative dimensionless ratios in the SU(3) gauge theory. 
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4.1 Technical comparison 

We can compare our results to those of [3]. We find^, 



srg = 11.6(1.1) 
trl= 7.1(1.4) 



this work 



15.8(3.2) 
7.5(2.8) 



Chen et al. 



(4.1) 
(4.2) 
(4.3) 



Comparing the matrix elements in units of tq saves us from having to settle on a value for 
ro m fm (ro ^ = 410MeV was used in [3]). We conclude that our results are in satisfactory 
agreement with those of [3] , and the statistical uncertainties have been reduced by a factor 
of at least two. 

We stress that while our continuum extrapolation uses data at lattice spacings (a/ro)^ 
in the range 0.010 to 0.035, the continuum extrapolation of [3] uses data in the range of 
spatial lattice spacings 0.044 to 0.22. The control of the continuum limit is thus qualita- 
tively different. One might worry that the support of the field strength operators used for 
the 'type IF discretization in [3] varies between 0.4fm and 0.9fm. On the other hand, the 
authors do find good agreement between certain observables which only become strictly 
equal in the continuum limit. 

We have used non-perturbative renormalization factors at every lattice spacing, while 
the procedure used in [3] for s and t implies that the continuum is approached asymptoti- 
cally with 0(^(0 ) as well as O(a^) corrections. 

Simulating at small lattice spacings comes at a heavy computational price for these 
observables: the signal-to-noise ratio for matrix elements of the gluonic dimension four 
operators decreases with the fourth power of the lattice spacing [2] . This explains why our 
statistical errors on the bare matrix elements are somewhat larger than in [3]. 

The use of a coarse spatial lattice allowed the authors of [3] to reach volumes signifi- 
cantly larger than ours in physical units. We have however checked for finite volume effects 
explicitly, and find no significant variation in the matrix elements. It is known [19] that the 
low-lying spectrum of glueballs exhibits remarkably small finite-size effects for L > 2.5ro 
(and for certain channels they remain small for even smaller box sizes). 

4.2 Glueball production rate in radiative J/ip decays 

First of all, we can compare our result for s with the QCD sum rule prediction [4] 



Using Go = 0.012GeV'^ (the 'gluon condensate') and r^^ = 410MeV this leads to sr§ 6.0, 
almost a factor two smaller than our result. Since we regard the gluon condensate as a 
phenomenological parameter, this disagreement does not surprise us too much. 

We may use our value for s to estimate the partial width of J/ ijj to radiatively decay into 
a scalar glueball. An approximate expression for this experimentally observable quantity 

^As compared to [3], our s contains an extra factor of ll/(47r)^ (coming from the beta- function). 




(4.4) 
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was derived in [4], 



Q 4 2 

rw^^Go7)^ ,,.„,.3,„grw7^^,^,) . (4.5) 



Here mc 1250MeV is the charm mass in the MS scheme at scale fx = rric, rrij/^ 2± 
3097MeV is the J/tp mass, a ~ ^ is the fine structure constant and T^J/ip — e'^e~) ~ 
5.55kcV. The ingredients that go into this formula arc the l/nic expansion, a dispersion 
relation for the charm-quark loop induced 77 gg transition and the assumption that the 
J/V' contribution dominates the spectral integral. This leads to the estimate 

Br( J/V' ^ G07) ~ 0.009. (4.6) 

The authors of [4] also suggest that a more accurate prediction is 

Here the assumption about the dominance of the J/V' is replaced by the assumption that 
the relative size of the J/^ and the 'continuum' contributions are identical in the scalar 
and the pseudoscalax channels. Using Eq. (5) of [20], we can trade the properties of rj' for 
the corresponding ones for 77 and use the SU (3)f result [20] 

{n\asF;,F;M = y yf /X (4.8) 
with fr, ITOMeV. Then we obtain 

Br(./* ^ Go.) . ^ (1^)' £^^I^^^ = 5.6(1.7) . B,(J/^ ,-,). 

(4.9) 

Here X = mrf/mj^,^ and x' = m^i /rriji^. Since Br(J/V' f]^) = 0.98(10) • 10^'^ [21], we 
obtain a somewhat lower result than (4.6). This is still a rather large branching ratio. 
For instance, the PDG [21] gives Br(J/V' ^ 7/o(1710)) 1.5(3) • 10"^ if one adds up the 
contributions of the 'yKK, jinr and 70^0; channels. The production rate of /o(1500) is 
even smaller. According to (4.9), if any of the /o(1370), /o(1500) or /o(1710) states had a 
glucball component close to unity, they would be produced more copiously than observed. 
Therefore, Eq. 4.9 suggests that the glueball component is quite strongly diluted among 
the three states. This conclusion is also reached by doing detailed parametrizations and 
fits to experimental data of the mixing pattern [22, 23]. 

4.3 Constraining thermal spectral functions 

Another application arises in the calculation of transport coefficients in the plasma of 
gluons. Indeed the local-current two-point function ^00 can be expressed in terms of the 
spectral function p(Ld,p,T), 

coshLj(2^ — r) 



^ /""^ . coshLj^^-r 

Aoo(r,p,r)=/ dwpiu;,p,T) . (4.10) 

Jo smn 
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Here T is the temperature. In the tensor channel the shear viscosity is then given by 

r,(r) = TT lim ^^^^^ . (4.11) 

Because ^oo is dominated by ultraviolet contributions that are temperature independent, 
in [24] we proposed to subtract from Aqq what that correlator would be if the spectral 
function was the same as at T = 0, before solving Eq. 4.10 for p{uj,p,T). Let us call that 
would-be correlator ^qo- The T = spectral function has a simple expression in terms of 
energy levels and matrix elements of the kind wc calculated in this paper, 

p{lo,p,T = 0) = L-'^J2^i^ - ^n)mOo\n)\\ (4.12) 

n 

Therefore, we find that the contributions of the first two terms in Eq. 4.12 to Aoo{t, 0, T) 
read (at Nt = l/{aT) = 8) 

Ioo(r, 0, 1.24Tc) = 1.25 + 5.2 + . . . (4.13) 
Ioo(t, 0, 1.65Ta) = 0.388 + 3.7 + . . . (4.14) 

These contributions are substantial since Ao(){t,0,T) = 8.05(32) and 8.73(33) respec- 
tively [24]. The first excited state contribution appears to be particularly large. Due 
to possible higher state contamination in our estimate of Fs* and Ft* (see section 3), it 
may be that the second terms in Eq. 4.13 or 4.14 effectively amounts to the contribution 
of more than one state. 

4.4 A summary of the low-energy parameters of SU(3) gauge theory 

We finish with a summary of our knowledge of low-energy dimensionless quantities in 
the SU(3) gauge theory. The results we use were all obtained with the Wilson plaquette 
action [8] and we extrapolate them to the continuum. Wc use the paramctrization of 
(ro/a)(/3) from [9]. If a is the string tension, Tc the deconfining temperature and the 
Lambda parameter in the MS scheme, we find in the continuum 



Msro 


= 3.958(47) 


XV(4 - 2) 


= 0.09, 


(4.15) 


MtVo 


= 5.878(77) 


XV(4 - 2) 


= 0.6, 


(4.16) 


srI 


= 11.6(1.1) 


XV(3 - 2) = 


= 0.5, 


(4.17) 


trI 


= 7.1(1.4) 


XV(3 - 2) = 


1.0, 


(4.18) 




= 1.1611(95) 


XV(5 - 2) 


1 = 0.1, 


(4.19) 


TcTo 


= 0.7463(64) 


xV(4-2) 


1 = 0.3, 


(4.20) 




= 0.60(5). 






(4.21) 



The data for aMs, qMt sa^ and ta'^ is from this work, plus the /3 = 6.1 glueball mass 
data from [10]. The string tension a^/a is taken from [17] {/3 > 5.8000) and from [10] 
(/? = 6.4). We took the critical values of P for given Nt = l/{aT) from [25] {Nt = 6, 8, 12) 
and from [26] {Nt = 5). The Lambda parameter value is the data of [27]. We do a 
continuum extrapolation linear in (a/ro)^ for the first five quantities. Wc extrapolate T^ro 
as a function of l/N^- The quantity Aj^ro is already given in the continuum in [27]. 
The overall level of accuracy achieved is remarkable. 



5. Conclusion 

In summary, we have computed scalar and tensor glueball matrix elements in SU(3) gauge 
theory with improved precision. A straightforward application to the production of scalar 
glueballs in J radiative decays suggests that none of the known scalar mesons can contain 
too large a glueball component. Finally we gave a summary of the current knowledge 
of seven non-perturbative observables of the SU(3) gauge theory. We hope that these 
quantities constitute a useful set for models and semi-analytical methods to calibrate on 
and compare predictions to. 
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